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ABSTRACT 

The amplitude of the Compton-Getting (CG) anisotropy contains the power- 
law index of the cosmic-ray energy spectrum. Based on this relation and using 
the Tibet air-shower array data, we measure the cosmic-ray spectral index to be 
— 3.03±0.55stat± < 0.62syst between 6 TeV and 40 TeV, consistent with —2.7 from 
direct energy spectrum measurements. Potentially, this CG anisotropy analysis 
can be utilized to confirm the astrophysical origin of the "knee" against models 
for non-standard hadronic interactions in the atmosphere. 

Subject headings: cosmic rays — solar system: general 



1. Introduction 

The cosmic-ray (CR) energy spectrum has been observed using ground-based air shower 
detectors and various instruments on board satellites and balloons. It features a power law 
dN/dE oc E~'^ in the energy range from 30 GeV to 100 EeV, where A^ denotes the number 
of cosmic rays at the top of the atmosphere, E the cosmic-ray particle energy. The spectral 
index 7 changes from 2.7 to 3.1 around the "knee" at E' ~ 4 PeV. Various models have been 
proposed for th e origin of this knee; most models suppose that it is of astrophysical origin 



( jHorandel 20041 ). while there are some models that attribute it to non-standard high-energy 



nucleon interactions in the Earth's atmosphere. 

For a power-law CR energy spectrum, the Compton-Getting (CG) anisotropy involves 
the spectral index 7. When a detector observing cosmic rays moves relative to the rest frame 
of the CR plasma in the Earth orbit around the Sun, a fractional CR intensity variation A/ 
can be observed in the solar time frame. For a presumed power-law CR energy spectrum. 
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this CG anisotropy (ICompton &: Getting 19351 : iGleeson fc Axford 19681 ) is expressed as 



^^ (7 + 2)-cos^, (i: 



< / > c 

where / denotes the CR intensity, 7 the power-law index of the cosmic-ray energy spectrum, 
V the orbital velocity, c the speed of light and 6 the angle between the arrival direction of 
CRs and the moving direction of the observer. 

At multi-TeV energie s, a clear evidence for t he CG anisotropy in the solar diurnal varia- 



tion was first reported in (ICutler &: Groom 19861 ). The observed variation was in reasonable 



agreement with a sinusoidal curve expected from the CG anisotropy with the maximum 
phase of the curve deviating from 6:00 h by +2 hours at 2a significance. This deviation 
was ascribed to meteorological effects on the underground muon intensity. The Tibet air- 
shower experiment performed the first study on energy dependence of the solar diurnal vari- 
ation and showed that i t was consistent with the CG anisotropy above multi-TeV energies 



(jAmenomori et al. 2004 : lAmenomori et al. 2006al ). 



Conventionally, the spectral index above ~100 TeV is obtained by indirect measurements 
with air-shower arrays, in which each CR particle energy is determined from its shower size. 
The relation between particle energy and shower size is deduced from air-shower simulations 
assuming hadronic interaction models as a reasonable extrapolation of the accelerator results 
to the higher energy range. The spectral index thus obtained is valid only when no violation 
of the standard hadronic interaction is expected in the relevant energy range. This fact 
leaves a room for an argument whether the knee is due to non-standard high-energy nucleon 
interactions in the atmosphere or not. 

On the other hand, the CG anisotropy does not depend on CR particle energy by Eq.([T]). 
We can then measure the power-law index of the all-particle energy spectrum in the way it 
is virtually uninfluenced by hadronic interaction models, because air-shower simulations are 
used only for the purpose of estimating the rough energy range where we are observing the 
CG anisotropy. Potentially, this method makes it possible to confirm that the bend of the 
energy spectrum at the knee is of astrophysical origin, not due to non-standard hadronic 
interactions in the atmosphere. 



2. Tibet Air-Shower Experiment 

The Tibet air-shower experiment has been operating successfully at 90.522° E, 30.102° N 
and 4300 m above sea level since 1990. Be ing upgraded several tirnes, the Tibet III array w e 



use here was completed in November 1999 (jAmenomori et al. 19921 : lAmenomori et al. 20021 ). 
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It consists of 533 scintillation counters of 0.5 m^ each. Each counter is viewed by a fast- 
timing (FT) photor aultiplier tube and is p laced on a 7.5 m square grid with an enclosed 



area of 22,050 m^ (jAmenomori et al. 20031 ). A 0.5 cm thick lead plate is put on top of 
each counter to improve fast-timing data by converting gamma rays into electron-positron 
pairs. An event trigger signal is issued when any fourfold coincidence takes place in the FT 
counters recording more than 0.6 particle, which results in the trigger rate of about 680 Hz 
at a few-TeV threshold energy. We estimate the energy of a primary CR particle by XIpft, 
which is the sum of the number of particles/m^ for each FT counter. 



3. Cosmic- Ray Data Analysis 

We collected 7.7 x 10^*^ events during 1319 live days from November 1999 to November 
2005. After some data selections (i.e., software trigger condition of any fourfold coincidence 
in the FT counters recording more than 0.8 particle in charge, zenith angle of the arrival 
direction less than 45°, and air shower core position located in the array), the events left 
were histogrammed in hourly bins in the local solar time according to their event time and 
incident direction. To suppress the seasonal change in the daily variation, we corrected the 
number of events for each month considering monthly variations of the live time. 



The daily and yearly event rates vary by ± 2% and ± 5% (jAmenomori et al. 20031 ) 



respectively, affected mostly by meteorological effects. To rem ove these temporal varia tions 



we adopted the following East— West subtraction procedure (JNagashima et al. 19891 ). We 
first obtain E{t) and W{t), the daily variation at solar time frame for each of East and 
West incident events, according to the geographical longitude of incident direction of each 
event. Subsequently dividing E(t) — W(t) by 26t, an hour-angle separation between the 
mean directions of E- and W-incident events, we finally reach D{t), the differential of the 
physical variation at solar time frame R{t). 

^ E(t) - W(t) ^ R(t + it)-R[t-U) ^ d 

^ ' 26t 26t dt ^ ' ^ ' 

The advantage of this method is that meteorological effects and possible detector biases 
which are expected to produce common variations for both E- and W-incident events largely 
cancel out. 
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4. Results and Discussion 

The CR events were divided into eiglit ^ pft bins in order to find out the energy region 
where the observed solar daily variation is free from effects other than the CG anisotropy. 
The fitting curve for the observed solar daily variation is expressed as 

/(A) = acos(^(A-0)), (3) 

where a denotes the amplitude of the CG anisotropy, A [hr] the local solar time and </> [hr] the 
phase where the sinusoidal curve reaches its maximum. Note that in a differential form, the 
phase is shifted earlier by 6 hours (1/4 cycle) from the corresponding actual daily variation 
and the amplitude is 7r/12 times that of the physical CG anisotropy. 

Figure [11(a) shows ^Pft dependence of the phase in Eq.(l3]). Although there exists 
deviation due to some contamination in the region of ^Pft < 50, the measured phase is 
consistent with the expected value of the CG anisotropy in ^ ppT > 50. Therefore, we 
regard the observed variation in ^ ppT > 50 as caused by the CG anisotropy. Figure [11(b) 
shows X^Pft dependence of the observed amplitude a. Since the CR energy spectrum in 
the multi-TeV region is expressed by a single power law, the amplitude should be consistent 
with a constant line if there is no contamination. This figure confirms that we do observe 
the CG anisotropy in the region of XIpft ^ 50. Thus, we report here on the spectral index 
using 9.7 X 10^ events with ^ ppT > 50. This region is roughly between 6 TeV and 40 TeV, 
corresponding to 20% and 80% of the median energy distribution, respectively. 

Figure [2]^a) shows the observed solar daily variation data fitted with a sinusoidal curve. 
The x^"fitting results are summarized in Table I. From the amplitude of the solar daily 
variation, we obtained 7 of CRs using the following relation, 

« = (7 + 2)^^F, (4) 

where a denotes the amplitude, v = 2.978 x 10^ [m/s] the average orbital velocity of the Earth 
and F = 0.827 an effective geometrical factor by which a decreases, calculated according to 
the latitude at Yangbajing site (see Appendix). The factor 7r/12 emerges, because we measure 
the anisotropy in the differential form Eq.([2]). As a result, we found a = {10. 9 ±1. 2 stat) x 10~^. 
In Eq.(j4]), we omitted a minor correction. Strictly speaking, the orbital velocity of the Earth 
is time dependent. Although it varies ±2%, we used the average v instead. It may become 
necessary to take this correction into account, when more data are accumulated. 

A seasonal change of the sidereal daily variation due to the galactic anisotropy could 
produce a spurious variation in the solar time frame. The differential variation in the local 
extended-sidereal time (367.2422 cycles/yr) is shown in Figure [2t^b). The insignificant varia- 
tion in the extended-sidereal time {0.7±1.2stat) x 10~^ ensures that contamination in the solar 
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daily variation due to the seasonal change of the sidereal daily variation is less than 12% of the 
CG anisotropy. Adding this as a systematic error, we found 7 = 3.08 ± 0.55stat± < 0-Q'^syst- 
This systematic error would become smaller in the future when more data are accumulated. 

As the detection efficiency of the Tibet air-shower array depends on CR nuclei in this 
energy range (6 — 40 TeV), the correction for it should be taken into account. An air- 
shower simulation employing direct obser vational data for primary cosmic rays revealed that 



0.05 should be added to the 7 above (jAmenomori et al. 2006bl ). This correction is not 



needed, however, above ~60 TeV where the detection efficiency becomes independent of 
nuclei, reaching ~100%. 

Finally, the spectral index 7 between 6 and 40 TeV turned out to be 3.03 ± O.bbstat^ < 
0.62syst, consistent with 2.7 from direct energy spectrum measurements. 

Using this CG anisotropy method, a future high-statistics experiment with a huge ef- 
fective area will be able to measure the power-law index of the all-particle energy spectrum 
at higher energies. The spectral index 7 measured by this method can be compared not 
only with those 7 values by direct measurements below ~100 TeV but also with those by 
indirect air-shower analyses above ~100 TeV. This would check the mutual consistency be- 
tween different ways to measure the CR energy spectrum. Following the same strategy, this 
CG anisotropy analysis can be applied to much higher energies to confirm the astrophysical 
origin of the knee. 

The collaborative experiment of the Tibet Air Shower Arrays has been performed under 
the auspices of the Ministry of Science and Technology of China and the Ministry of Foreign 
Affairs of Japan. This work was supported in part by Grants-in-Aid for Scientific Research 
on Priority Areas (712) (MEXT), by the Japan Society for the Promotion of Science, by the 
National Natural Science Foundation of China, and by the Chinese Academy of Sciences. 



A. APPENDIX 

Coefficient F in Eq.(jl]) is a correction factor for the latitude at Yangbajing site. To 
calculate F, cos 6 in Eq.([T]) {6: the angle between the arrival direction of CRs and the 
moving direction of the observer) appears as 

cos^ = cos L cos ip sm{ip + ip) — sin T sin L sirup 

— cosTcosLsin(y9cos(v9 + ^) (Al) 

with ij = 27rt/24, ip = 27rt/24/365.2422, 
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where L = 30.102° denotes the latitude at the Tibet air shower array, T = 23.44° the tilt 
angle of the Earth's spin axis relative to the orbital axis, t [hr] the lapse time from the 
beginning of a year, ip the solar time and if a variable which advances from to 27r in a 
year. Suppose the Earth's spin axis were perpendicular to the ecliptic plane (T = 0°) for 
simplicity, Eq. lJAip would reduce to cos^ = cos L sin ^. This means the amplitude of CG 
anisotropy decreases by a factor F = cosL = 0.865. In reality, however, we need to take 
into account the tilt of the Earth's spin axis according to Eq. OAip . We found F = 0.827, 
approximately 4% smaller than 0.865, by numerically integrating Eq. OAip so that cos 6 can 
be expressed in terms of the solar time ip only. 
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Table 1: The x^-fitting results of the differential variation in the local solar time (a) and in 
the local extended-sidereal time (b), assuming the sinusoidal curve Eq.([n]). The error bars 
are statistical only. 



a (x 10-3%) [hr] xVd.o.f. 

10.9 ±1.2 -0.14 ±0.41 8.03/22 



(b) 0.7 ±1.2 -4.0 ±6.9 21.9/22 



-10- 



2 

I ' 

; » 

ft 

-2 
-3 













1 


; 


+ * t 








[] 




1 










1 


(a) 




: 



10 



100 
ZpFT 



1000 



0.025 

0.02 

b 0.015 

0.01 






0.005 




[] 






T 






+ 1 
" 11 




1 














1 


(b) 







10 



100 
EppT 



1000 



Fig. 1. — (a) : XIpft dependence of the phase 0, (b) : XIpft dependence of the amphtude 
a. Both are evaluated in a differential form. The solid lines show the expected values from 
the CG anisotropy, assuming 7 = —2.7. The error bars are statistical only. 
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Fig. 2. — (a) : The differential variation of relative CR intensity in the local solar time -D(t), 
(b) : the differential variation of relative CR intensity in the local extended-sidereal time. 
The solid lines are the sinusoidal curves best fitted to the data. The error bars are statistical 
only. 



